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Abstract. This paper investigates small-scale (500 AU) structures of dense gas and dust around the low-mass pro- 
tostellar binary NGC 1333-IRAS2 using millimeter-wavelength aperture-synthesis observations from the Owens 
Valley and Berkeley-IUinois-Maryland- Association interferometers. The detected A = 3 mm continuum emis- 
sion from cold dust is consistent with models of the envelope around IRAS2A, based on previously reported 
submillimeter-continuum images, down to the 3", or 500 AU, resolution of the interferometer data. Our data 
constrain the contribution of an unresolved point source to 22 mjy. The importance of different parameters, such 
as the size of an inner cavity and impact of the interstellar radiation field, is tested. Within the accuracy of the 
parameters describing the envelope model, the point source flux has an uncertainty by up to 25%. We interpret this 
point source as a cold disk of mass > 0.3 Mq. The same envelope model also reproduces aperture-synthesis line 
observations of the optically thin isotopic species C^'^S and H^^CO"^. The more optically thick main isotope lines 
show a variety of components in the protostellar environment: N2H[^ is closely correlated with dust concentrations 
as seen at submillimeter wavelengths and is particularly strong toward the starless core IRAS2C. We hypothesize 
that N2H"'" is destroyed through reactions with CO that is released from icy grains near the protostellar sources 
IRAS2A and B. CS, HCO^, and HCN have complex line shapes apparently affected by both outflow and infall. In 
addition to the east-west jet seen in SiO and CO originating from IRAS2A, a north-south velocity gradient near 
this source indicates a second, perpendicular outflow. This suggests the presence of a binary companion within 
0'.'3 (65 AU) from IRAS2A as driving source of this outflow. Alternative explanations of the velocity gradient, 
such as rotation in a circumstellar envelope or a single, wide-angle (90°) outflow are less likely. 
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1. Introduction 

Our understanding of the cloud cores that form stars 
has benefited significantly from the advent over the last 
years of (sub)millimeter-continuum bolometer cameras. 
Sensitive, spatially resolved measurements have allowed 
quantitative testing of models of starl ess /pre-stellar cores 
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of the matter distribution during star formation, they 
also can serve as 'baselines' for interpreting higher res- 
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olution observations obtained with millimeter interferom- 
etry. Such data address the presence and properties of 
circumstellar disks during the e arly embedded pha se (e.g. 
iHogerheiide et a,lJll99Sl llflQfll: iLoonev et al.' 2000). This 
Paper presents millimeter aperture-synthesis observations 
of continuum and line emission of the young protobi- 
nary system NGC 1333-IRAS2, and uses modeling re- 
sults bas ed on single-dish submillimeter continuum imag- 
ing from lJorgensen et al.l (2002; Paper I hereafter) to in- 
terpret the data in terms of a collapsing envelope, a disk, 
and (multiple) outflows on 500 AU scales. 

The dee ply embedded ('class '; iLadal Il987l: 
lAndre et all EsSS) young stellar system NGC 1333- 
IRAS2 (IRAS 03258-f3104; hereafter IRAS2) has been 
the subject of several detailed studies. It is located in 
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the NGC 1333 molecular cloud, well known for harboring 
several class and I obiects. and was first identified 
from IRAS data by Ijennings et al.l lll987|). Qu oted dis- 
tances t o NGC 1333 range fro m 220 PC l)Cernia.l99 0) to 
350 pc i Herbig fc JonejFlQSSl) : here we adopt 220 pc in 
rdance with lPapCTl At this distance the bolometric 



Table 1. The parameters for IRAS2 from lPaoeTl . 



accorc 



luminosit y of IRAS2 is 16 Lf^. Submillimeter-continuum 
imaging llSandell fc l200ll and Fig. m below) and 

high-resolution mi llimeter interferometry IjBlakd llQQf* 
iLoonev et alJl200fH) have shown that IRAS2 consists of at 
least three components: two young stellar sources 2A and, 
30" to the south-east, 2B; and one starless condensation 
2C, 30" north-west of 2A. The s ources 2A and 2B are 
also detected at cm wavelengths (jRodriguez et al.lll999t 
iReipurth et"al]l2002() . 

Maps of CO emission of the IRAS2 region show 

two outflow s, dir ected north-south and east-west 

(Liscau et all Il9 88": 'SandcU ct al.' 'l994; 'K nee fc SandelJ 
I2OOO.: .Engargiola fc Plambeck .199^) . Both flows ap- 
pear to originate to within a few arcseconds from 2A 
I Engargiola fc Plamb eck Il999() . indicating this source is 
a binary itself although it has not been resolved so-far. 
The different dynamical time scales of both flows suggests 
different ev olutionary stage s for t he binary members, 
which lead iKnee fc Sandelll l|2000l) to instead propose 
2C (30" from 2A) as driving source of the north-south 
flow. It is unclear how well dynamic time scales can be 
estimated for outflows that propagate through dense and 
inhomogeneous clouds such as NGC 1333. Single-dish CS 
and HCO"'" maps als o show contributions by the o utflow, 
especially for CS ijWard-Thompson fc Buckie^ EoOLl . 
The north-south outflow may connect to an observed 
gradient in centroid velocities near 2A, but the authors 
cannot rule out rotation in an envelope perpendicular to 
th e east- we st flow. 

IPaoer j determined the physical properties of the 
IRAS2 envelope using one-dimensional radiative trans- 
fer modeling of Suhmillimeter Common User Bolometer 
Array (SCUBA) maps and the long-wavelength spec- 
tral energy distribution (SED). Assuming a single radial 
power-law density distribution, p cx r"^, an index p = 1.8 
and a mass of 1.7 Mq within 12,000 AU was found (see 
Tableland Fig.^. Montc-Carlo modeling of the molecu- 
lar excitation and line formation of C^*0 and C^^O obser- 
vations yield a CO abundance of 2.6 x 10^^ with respect 
to H2, a fac tor 4-10 lower than w hat is found in local dark 

Il982t l 



clouds fe.g. lFrerking et alJII 982HT>a,cv et a]]l^994^ . 

This Paper presents A = 3 mm interferometric obser- 
vations of IRAS2 in a range of molecular emission lines 
probing dense gas and continuum emission tracing cold 
dust. It builds on the modeling of Paper I by using it as a 
framework to interpret the small-scale structure revealed 
by the aperture-synthesis data. Section El describes the 
observations and reduction methods. Section |3| analyzes 
the continuum emission, and compares it to the previ- 
ously derived models. Section 0] presents the molecular- 
line maps and discusses the physical and chemical prop- 
erties of the gas in the proximity of IRAS2. Section re- 



Distance, d 


220 pc 


ibol 


16 Lq 


Tbol 


50 K 


Envelope parameters: 


Inner radius (T = 250 K), 


23.4 AU 


Outer radius, -Riok" 


1.2x10" AU 


Density at 1000 AU, n(H2) 


1.5x10*^ cm"^ 


Slope of density distribution, p 


1.8 


Mass, Ml ok" 


1.7 Mq 


CO abundance, [CO/H2] 


2.6x10"^^ 



Notes: "The outer boundary is not well constrained, but taken 
to be the point where the temperature in the envelope has 
dropped to 10 K. The mass refers to the envelope mass within 
this radius. 



ports a pronounced north-south velocity gradient around 
IRAS2A and explores rotation or outflow as possible ex- 
planations. Scction|Slconcludcs the Paper by summari zing 
the m ain findings. A companion paper ( J0rgcnscn ct alJ 
l2003al) presents a detailed study of the bow shock at the 
tip of the east-west jet from IRAS2 based on single-dish 
and interferometric (sub)millimeter observations. 

2. Observations 

2.1. Interferometer data 

IRAS2 (a(2000) = 03'^28'"56':29; (5(2000) = 3in4'33"93) 
was observed with the Millimeter Array of the Owens 
Valley Radio Observatory (OVRO)^ between October 5, 
1994 and January 1, 1995 in the six-antenna L- and H- 
configurations. Tracks were obtained in two frequency set- 
tings at 86 and 97 GHz, and each track observed alter- 
natingly two fields: the source positions discussed in this 
Pap er and the bow shock at the end of the eastern jet 
( Jorgensen et al. 2003al. The observed tracks cover pro- 
jected baselines of 3.1-70 kA at 86 GHz. The observed lines 
are listed in Table [S] and were recorded in spectral bands 
with widths of 32 MHz (- 100 km s^^). Y{^^CO+ 1 - 
and CS 2—1 were observed in 128 spectral channels and 
the remaining lines in 64 spectral channels. The complex 
gain variations were calibrated by observing the nearby 
quasars 02344-285 and 3C84 approximately every 20 min- 
utes. Fluxes were calibrated by observations of Uranus 
and Neptune. Calibration and flagging of visbilities with 
clearly deviating amplitudes and/o r phases was perform ed 
with the MMA reduction package ijScoville et al.lll99^ . 

The millimeter interferometer of the Berkeley-Illinois- 
Maryland Association (BIMA)^ observed IRAS2 on 
November 4-5, 2000, and January 20-21, February 20, and 
June 5-6, 2001. The array B-, C-, and D-configurations 

^ The Owens Valley Millimeter Array is operated by the 
California Institute of Technology under funding from the US 
National Science Foundation. 

^ The BIMA array is operated by the Universities of 
California (Berkeley), Illinois, and Maryland, with support 
from the National Science Foundation. 
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Fig. 1. a) and b) SCUBA maps of IRAS2 at 450 and 850 ^m, respectively, centered on IRAS2A. c) Observed SED 
(symbols) and model fit (solid cm've). d) and e) Brightness profiles at 450 and 850 /im (symbols) and the model 
fit (solid curve). The dashed curves show the beam profiles, f) Density (dashed line) and temperature (solid line) 
distributions of best-fit model. See Paper I for details. 



provided projected baselines of 1.7-68 kA. The lines of 
HCO+ J=l-0, HCN 1-0, N2H+ 1-0, and C^^S 2-1 were 
recorded in 256-channel spectral bands with a total width 
of6.25MHz (~ 20 kms~^).The complex gain of the inter- 
ferometer was calibrated by observing the bright quasars 
3C84 (4.2 Jy) and 0237-^288 (2.3 Jy) approximately ev- 
ery 20 minutes. The absolute flux scale was bootstrapped 
from observations of Uranus. The rms noise levels are 
0.14 Jy beam"-'^ in the 24 kHz channels, with a synthe- 
sized beam size of 8.2" x 7.5" FWHM. The data were cal- 
ibrated with routi nes from the MIRIAD software package 
i Sault et al.1ll995|) . 



In the reduction, data points with clearly deviating 
phases or amplitudes were flagged. The maps were cleaned 
down to 3 times the rms noise using the MIRIAD 'clean' 
routine. The strong continuum of the two central point 
sources allowed self-calibration, which was applied and 
used to correct the spectral line data. The naturally 
weighted continuum observations typically had rms noise 
better than lxlO~^ Jy beam~^ with half power beam 
widths (HPBW) of « 3" for the OVRO observations and 
« 8" for the BIMA data (see Table E}. Table lists the 
details of the line observations. 



Table 2. Line data of IRAS2 discussed in this Paper. 



Molecule 


Line 


Rest freq. 


Observed with 


CH3OH 


2i - li 


97.5828 


OVRO 


CS 


2-1 


97.9810 


OSO, OVRO 




3-2 


146.9690 


IRAM 30m 




5-4 


244.9356 


IRAM 30m, JCMT'' 




7-6 


342.8830 


JCMT'' 




2-1 


96.4129 


IRAM 30m, BIMA 




5-4 


241.0161 


JCMT 


HCN 


1-0" 


88.6318 


OSO, BIMA 




1-0 


86.7543 


OSO, OVRO 


HCO+ 


1-0 


89.1885 


OSO, BIMA 


N2H+ 


1-0° 


93.1737 


OSO, BIMA 


SiO 


2-1 


86.8470 


OVRO 


so 


22 - li 


86.0940 


OVRO 


S02 


73,5 — 82,6 


97.7023 


OVRO 



Notes: "Hyperfine splitting observed in one setting. ''Archival 
data. 
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Table 3. Results of fits to the visibilities. 



OVRO BIMA 


RMS (Jy beam"') 
Beam 


1.0x10"^ 0.9x10"-' 
3.2" X 2.8" 8.2" X 7.5" 


IRAS 2 A 


Kot (Jy) 
X-ofTset (") 
Y-offset (") 


0.035 0.040 
-9.17 -8.92 
3.67 3.10 


IRAS2B 


Ftot (Jy) 
X-offset (") 
Y-offset (") 


0.012 0.014 
13.66 13.78 
-17.51 -18.37 



Notes: IRAS2A is marginally resolved, whereas IRAS2B is un- 
resolved. 



2.2. Single-dish data 

Using the Onsala 20 m telescope (OSO)-^ a number of 
molecules were observed towards IRAS2A in March 2002. 
These are listed in T a ble El We al so include line data 
from |j(3rgensen et all l|200a l2003hl) obtained with the 
IRAM 30m4 and James Clerk Maxwell (JCMT)^ tele- 
scopes, and spectra taken from the JCMT archive^. The 
lines were converted to the main beam antenna temper- 
ature scales using the appropriate efRciencies, and low- 
order polynomial baselines were fitted and subtracted. 



3. The continuum emission 

The 3 mm continuum images clearly show the two compo- 
nents 2A at (-9", +3") and 2B at (-f 14", -18") (Fig.©. 
Table|311ists the results of fits of two circular Gaus sians to 
the visibility data. Consistent with iLoonev et al\ (|200G), 
2A is the stronger of the two. Differences in detected 
fluxes between the OVRO and BIMA data sets indicate 
that the emission is extended, and varying amounts are 
picked up by the respective (u, v) coverages of the arrays. 
Emission from the third source 2C, north-west of 2A, is 
not detected, supporting the suggestion that it has not yet 
formed a star and lacks a strong central concentration. 



^ The Onsala 20 m telescope is operated by the Swedish 
National Facility for Radio Astronomy, Onsala Space 
Observatory at Chalmers University of Technology. 

The IRAM 30 m telescope is operated by the Institut de 
Radio Astronomic Millimetrique, which is supported by the 
Centre National de Recherche Scientifique (France), the Max 
Planck Gesellschaft (Germany) and the Instituto Geografico 
Nacional (Spain). 

^ The JCMT is operated by the Joint Astronomy Centre in 
Hilo, Hawaii on behalf of the parent organizations PPARC in 
the United Kingdom, the National Research Council of Canada 
and The Netherlands Organization for Scientific Research 

® The JCMT archive at the Canadian Astronomy Data 
Centre is operated by the Dominion Astrophysical Observatory 
for the National Research Council of Canada's Herzberg 
Institute of Astrophysics. 
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RA offset ["] RA offset ["] 

Fig. 2. Maps of continuum emission at 86-89 GHz from 
OVRO (a) and BIMA (b). Offsets are with respect to the 
pointing center of a(2000) = 03^28'°56!29 and (5(2000) = 
31°14'33"93. Contours are shown at 3a, 6a, 12a, 24a and 
48cr, where a is the RMS noise level of Table The filled 
ellipses in the lower left corner of the panels indicate the 
synthesized beam sizes; the large circles show the 50% 
sensitivity levels of the primary beams. 

3.1. A model for the continuum emission 

The different detected fluxes from OVRO and BIMA in 
Table |3 and comparison of the interferometer images of 
Fig. 121 and the SCUBA images of Fig. ^ clearly show that 
the arrays have resolved out significant amounts of ex- 
tended emission because of their limited {u, v) coverage. 
The envelope model (density, tem perature, dust emissiv- 
ity as function of wavelength) by IPaner j predicts sky- 
brightness distributions at 3 mm, and fluxes in the inter- 
ferometer beams after sampling at the actual {u, v) po- 
sitions and subtracting the contribution from 2B. This 
latter subtraction of the Gaussian fit to 2B only affects 
the results minimally, indicating that 2B is well separated 
from and much weaker than 2A. Fig. |3 compares the pre- 
dicted flux as function of projected baseline length with 
the data. In addition to the model envelopes, we have in- 
cluded as free parameter the flux of an unresolved point 
source (< 3"). Because a point source contributes equally 
on all baselines, this addition corresponds to a vertical 
offset of the model curve in the plots. Such offsets are ap- 
parent in both OVRO and BIMA data, and a point source 
flux of 22 mjy at 3 mm provides an adequate fit to the 
data when added to the model envelopes. 

The infer r ed p oint source flux is consistent with 
iLoonev et al.l l(2003l) . who flnd a 20 mJy source at 2.7 mm 
associated with IRAS2A. For reasonable assumptions 
about the spectral index of the point source (q;=2-4 if 
thermal) this component does not contribute signiflcantly 
to the flux in the much larger SCUBA beams, and there- 
fore does not invalidate the SCUBA-based envelope mod- 
els. The thermal nature of the point source is supported by 
detection of 2A at a flux of 0. 22 mJy at 3.6 cm and a res- 
olution of 0.3" with the VLA l)Reipurth et al.ll2002|) . This 
yields a spectral index of 1.9 between 3.6 cm and 3.3 mm, 
consistent with optically th ick thermal emiss ion. A similar 
conclusion was reached bv iRodriguez et sl\ (jl99& 'l based 
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Fig. 3. Visibility amplitudes of the observed continuum 
emission from the BIMA (upper panel) and OVRO (lower 
panel) observations as a function of projected baseline 
length in kA centered at the position of IRAS2A. The 
data are plotted as filled symbols with la error bars, the 
continuous lin es indicat e the predictions from the contin- 
uum model of IPaoer^ with the same {u,v) sampling as 
the observations - respectively with (solid) and without 
(dashed) an unresolved compact source of 22 mJy added 
to the model. The dotted histogram indicates the zero- 
expectation level: the expected amplitude signal due to 
noise alone in the absence of source emission. 



on the spectral index from VLA observations of IRAS2A 
at 3.6 and 6 cm. 

Assuming that the point source emission is optically 
thin and thermal, the inferred flux of 22 mJy corresponds 
to a dust mass of 3.3 x 10"'^ Mq if we adopt an average 
temperature of 30 K and a emissivity per unit (dust) mass 
at 3.5 mm o f k = 0.24 cm^ frq r n extr apolation of the 
opacities bv lOssenkonf fc Hennind ljl994l) for grains with 
th in ice ni antles as was assumed in the envelope models 
in IPaper"^ With a standard gas-to-dust ratio of 100, the 
total mass is 0.33 Mq. If the emission is optically thick as 
the spectral index indicates, this is in fact a lower limit to 
the mass. The favored explanation for this compact mass 
distribution is a circumstellar disk. 



3.2. Parameter dependency of the continuum model 

To test the validity of the envelope model a number of 
parameters were varied within the constraints set by the 
modeling of the SCUBA observations (Fig. The un- 
certainty in the power-law index from the SCUBA model 
of Paper I is ±0.2. Over this range of density slopes we 
find central point source fluxes of 22 ± 4 mJy, with a clear 
degeneracy between the slope of the density profile and 
the fl trx of the central p oint source. This is similar to 
what iHarvev et al.1 l|2003fl find in a detailed analysis of 
high-resolution millimeter continuum observations of the 
class object B335. Both BIMA and OVRO data sets 
are fitted well within the uncertainties using the density 
profile slope from the SCUBA data, although the actual 
best fit model to the OVRO data has a slightly steeper 
density slope {p — 1.9) and a lower point source flux. The 
interferometry data cannot constrain the slope of the den- 
sity profile further than its uncertainty from the SCUBA 
model. So although the data can be fitted with a single 
power-law density envelope from the scales probed by the 
SCUBA observations down to the scales probed by the 
interferometry observations, a steepening or flattening of 
the density profile at small scales cannot be ruled out. 

iHarvev et alJ (j2003} find that the uncertainty in the 
central point source flux dominates over uncertainties in 
other model parameters such as external heating by the 
interstellar radiation field (ISRF), wavelength dependence 
of the dust emissivity, outer radius of the envelope, and de- 
viations from spherical symmetry (e.g., an evacuated out- 
flow cavity) . Because our interferometer data only sample 
the inner regions, they are not sensitive to variations in 
the outer radius or inclusion of heating by the ISRF. The 
latter hardly affects the temperature structure because the 
source is relatively luminous and dominates the heating as 
is seen in Fig. |5| 

In simila r stud ies to that pre s ented in IPaoer 1 
IShirlev et all l)2002() and lYoung et all l)2003|) modeled the 
SEDs and brightness profiles from SCUBA observations 
of protostellar sources using ID radiative transfer, assum- 
ing power-law density profiles and solving for the tem- 
perature structure. Two differences exist, howev er, in the 
approaches t aken in these two papers and our IPaoer j 
IShirlev et alJ and lYoung et alJ included contributions to 
the heating of the envelope by the external interstellar 
radiation field and adopted outer envelope radii signifi- 
ca ntly larg er than those set by the 10 K boundary used 
in^P aperl F or the sources common to the two samples, 
yShirlev et alJ found on average steeper de nsity profiles 
than ours for the class objects whereas lYoung et alJ 
found similar dens ity profiles to ours for the class I ob- 
jects. lYoung et al. suggested that the disagreement for 
the class objects and agreement for the class I objects 
was due to a combination of neglect of the ISR F and an 
underestimate of the sizes of the envelopes in IPaoer ^ 
while inclusion of the ISRF will indeed tend to steepen 
the derived density profile, an overestimate of the outer 
radius (by factors of 2 or more) will tend to flatten the 
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Fig. 4. Visibility amplitudes of the observed BIMA continuum emission as in Fig. 13 compared to various input models 
centered at the position of IRAS2A. Upper panels: a) models with changing steepness of the density profile, b) test 
of different values of the outer radius and inclusion of the interstellar radiation field. Models with an outer radius 3 
times la rger than t he model from paper I (i.e. 36000 AU) are shown. Lower panels: c) fit to the inside-out collapse 
model of lShul l)l977j) with parameters constrained independently by molecular line observations and SCUBA continuum 
observations, d) models with changing size of the inner radius. 



derived profile. As illustrated above, however, these pa- 
rameters have negligible impact on the IRAS2 envelope 
structure. It is therefore interesting to note the agree- 
ment in slope between the interferometer and SCUBA 
continuum observation s, in contrast to the discussion of 
R335 bvlHarvev et al J E)B\. _Comparing to the results of 
IShirlev et all l|2002l) . iHarvev et alJ found a slightly flatter 
density profile when modeling the interferometer observa- 
tions. While uncertainty in the outer radius and ISRF may 
lead to only small departures for the interferometry data, 
it can lead to systematic changes in the slope of derived 
power-law density profile from the SCUBA observations 
of ~ 0.2. This could explain the differences between the 
density profiles from the interferometry and SCUBA data 
for B335. 

Our envelope model is entirely based on SCUBA data, 
and the interferometer fluxes serve only to constrain any 
point source flux. The robustness of that constraint de- 
pends on the assumption that the envelope model can 



be extrapolated down to scales much smaller than the 
SCUBA resolution (4"= 900 AU). In ' Paper ll the inner 
radius is fixed at a temperature of 250 K which occurs 
at i? = 22 AU, but it was argued that this is not deter- 
mined by the data. The size of any inner cavity is expected 
to affect the interferometer data since these sample small 
scales: a larger adopted cavity would result in a higher in- 
ferred point source flux to compensate for the reduced 
small-scale emission. Fig. plots this 'required' point 
source flux against inner cavity size. The point source flux 
increases from 22 mJy for cavities < 25 AU to « 27 mJy 
for cavities > 200 AU (!")■ The interferometer data would 
resolve cavities larger than this and a point source could 
no longer compensate for the removed emission. 

Turning this reasoning around, the inferred point 
source could be due to an increase in envelope density on 
small scales, as opposed to a circumstellar disk in an enve- 
lope cavity. Assuming a temperature of 150 K appropriate 
for the envelope scales unresolved by the interferometers 
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Fig. 5. Temperature profile in the outermost region of the 
envelope without (solid line) and with (dashed line) con- 
tributions from the interstellar radiation field. The dotted 
line indicates the temper ature of 1 K corresponding to 
the envelope outer radius ijPaper if) . 
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tions of IChiang Sz GoldreichI l)l997|) and iDullemond et alJ 

(200]]), and the envelope's inner cavity has been increased 
to 200 AU in radius so that it encompasses the disk. As 
a result, the temperature at the inner edge of the enve- 
lope drops from the original 250 K (at 22 AU) to 75 K 
(at 200 AU). The radiative transfer code DUSTY pro- 
duces the envelope's temperature distribution and emer- 
gent SED using the star-fdisk spectruin as heating in- 
put, similar to the calculations of 'Paoer J for the star-only 
spectrum. The comparison in Fig. [71 shows that the SEDs 
between 60 /im and 1.3 mm are unchanged. Our derived 
envelope parameters are therefore unaffected by the exact 
form of the input spectrum. The departures grow larger at 
the shorter wavelengths (2-20 /zm) and may be observable 
with, e.g., SIRTF. It is not surprising that the SEDs are 
most different at these wavelengths. Flared disk models 
such as those of IChiang fc GoldreichI l|l997ll are specifi- 
cally invoked to explain so-called 'flat-spectrum' sources. 
Their superheated surface layers 'flatten' the SED of these 
star-|-disk systems by boosting the 2-20 /im emission. It 
is not obvious that such a description of the disk is valid 
for early, deeply embedded objects, such as IRAS2A. Still, 
the important point here is that the influence of the disk 
on the observed SED is likely to be negligible at the wave- 
lengths where the envelope model is constrained. 
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Fig. 6. Derived point source flux plotted against size of 
the inner envelope cavity. 

instead of the 30 K assumed for the disk, one derives a 
mass of 0.06 AIq. For comparison the mass of the single 
power-law density model within 150 AU is only 0.008 Mq. 
So to explain the detected flux an increase in density by 
almost an order of magnitude is needed, which seems un- 
likely. 

The model from lPaper j assumes that the envelope is 
heated by a stellar blackbody of 5000 K at the center. If, 
as is argued above, the star is surrounded by a disk that 
reprocesses a significant fraction of the stellar light, the 
input spectrum shifts to longer wavelengths. To investi- 
gate the effect on the envelope's temperature structure, 
Fig. 13 compares the SEDs of the original model and a 
model where the central star is surrounded by a 200 AU 
outer radius, 0.33 Mq disk. The disk follows the descrip- 



3.3. A collapse model for the continuum emission 



As demonst r ated by iHogerheiide fc Sandelll l)2000|) . 
IShirlev et all l(2002f) . and ISchoier et alJ l(20o1 models 
other than a density power law can also fit continuum ob- 
seryations, in particular the inside-out collapse model of 
These authors conclude that a collapse model 



can provide an equally good fit as p ower-law models, with 
the caveat by IShirlev et alJ l)2002|) that collapse models 
only fit their class objects for sufficiently low ages where 
this model is well approximated by a single p ower law on 
the scales resolved by SCUBA. ISchoier et al.i (^02) and 
IShirlev et al. ( 2002) find that continuum and line data sets 
sometimes give discrepant collapse model fits to the same 
sources, with line data favoring higher ages than contin- 
uum data. 

Fitting the'Sb] (Il977l) inside-out collapse model to the 
SCUBA data for IRAS2 gives best fit values of a = 0.3 km 
s~"^ and t — 1.7 X 10"* years (see Fig.lSJ with the quality of 
the fits essentially identical to those of the single power- 
law models. This collapse model also fits the BIMA and 
OVRO data if a point source of 25 mJy is introduced. 

The integrated CS and C^'^S line intensities were fit- 
ted independently with the collapse rn odel using detailed 
radiative transfer as in ISchoier et al.l (|2002) for a con- 
stant fractional abundance with radius. The uncertainties 
in the line fluxes are assumed to be 20% and for each 
model the x^-estimator is used to pick out the best model 
and estimate confidence levels for the derived parameters. 
Interestingly, the fits to the CS and C^^'S fines (Fig.EJ give 
identical parameters to those derived from the dust mod- 
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Fig. 7. Changes of the emerging SED due to mclusion of 
a 200 AU outer, 0.3 Mq disk: in the upper panel the SEDs 
from the star + envelope (solid line) and star + envelope 
+ disk models (dashed line) are compared. In both models 
the central star is represented by a 5000 K blackbody. In 
the lower panel the ratio between the models are shown - 
the typical 20% error-level in the flux calibration is illus- 
trated by the solid rectangle. 
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Fig. 8. Fits to the SCUBA observations with a inside- 
out collapse model with an isothermal sound speed, a, of 
0.3 km s~^ and an age of 1.7 x 10'' years. 
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Fig. 9. Constraints on the inside-out collapse model de- 
rived from the CS and C'^^S line intensities, assuming CS 
and C^'^S abundances of Ix 10^^ and Ix 10^^° respectively 
found in the single power-law density model. The inferred 
age and sound speed agree well with those derived from 
the continuum data. 



eling, in c ontrast with the other sources fe.g. lShirlev et all 
2002; SchiereEaOHool. 

The identical fits to the line intensities and continuum 
observations and success of both collapse and power-law 
density models illustrates the low age inferred for IRAS2: 
the collapse expansion radius in the inside-out collapse 
model is for the envelope parameters located at « 4" 
and therefore not directly probed by the SCUBA contin- 
uum maps. Likewise the CS lines predominantly probe the 
outer regions of the envelopes where the density profile in 
the collapse model for the low age of IRAS2 essentially 
is a power-law. This also explains why a slightly higher 
point source fiux is obtained with the collapse model than 
the power-law density model: inside the collapse expansion 
radius the density profile fiattens in the collapse model, 
lowering the mass and thereby the fiux towards the un- 
resolved center of the envelope. This is compensated by 
increasing the point source flux when modeling the inter- 
ferometer observations. 

In summary, the interferometer continuum data are 
well described by the same 1.7 Mq envelope models that 
fit the SCUBA data. Power-law descriptions for the den- 
sity and inside-out collapse model fit the data equally 
well. They indicate the presence of a Fa^m ~ 22 mJy 
point source, presumably a > 0.33 Mq circumstellar disk. 
Uncertainties associated with the envelope model are re- 
fiected in the accuracy of the point source fiux, which may 
vary by up to 25% from the quoted value. The next section 
describes the line emission in this context. 
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4. Line emission 

4.1. Morphology 

Fig . 1 1 Ul iTl shows the integrated intensity maps of all lines 
detected with BIMA (HCN, HCO+, N2H+, C^^S) and 
OVRO (CS, H"CO+, SO, and CH3OH). Velocity centroid 
images are shown for CS, HCN and HCO"*" in Fig. [T^ 
Emission of SiO and SO2 was not detected toward the 
source position. The images from BIMA show more ex- 
tended structure than those from OVRO because of the 
different {u,v) coverage of the two arrays. All detected 
lines have a peak near the object 2A, and most show a 
peak near 2B. The non-detections of SO and CH3OH near 
2B are likely due to limited sensitivity given the low signal- 
to-noise of these lines toward 2A. Only the non-detection 
of N2H+ toward 2B is highly significant: emission in this 
line appears to avoid 2B. The extended emission picked 
up by BIMA shows three components. A roughly north- 
south ridge seen in HCN, HCO"*", and N2H+; emission 
along the east-west outflow in HCO"*" and, at the tip of 
the eastern jet at the edge of the image, in HCN; and an 
extended peak in N2H"'" around the continuum position 
2C. Interestingly, N2H+ also seems to avoid the east-west 
outflow and appears almost anti-correlated with HCO+. 

The intensity ratio of 1.5:2.7:4.7 of the N2H+ hyper- 
fine lines suggest that the emission is optically thin and 
close to LTE, where the ratio would be 1:3:5. Relative to 
the 450 /zm emission from SCUBA that traces cool dust, 
the N2H"'" emission is strongest around 2C, lower around 
2A, and absent toward 2B as illustrated in Fig . [13 1 In a 
study of dark cloud cores, iBers-in et al J ll2nnili2nni find 
that N2H"'" has a large abundance deep inside the clouds 
and a lower abundance in the exterior regions. This trend 
is opposite to tha t of CO, whic h is often highly depleted 
deep inside cores. iBergin et al.l argue that N2H+ is effec- 
tively destroyed through reactions with CO, and therefore 
is only present at high abundance where CO is depleted. 
This scenario can also explain the relative distribution of 
N2H+ in 2C, 2A, and 2B: in the starless core 2C tem- 
peratures are low and CO is highly depleted, resulting in 
strong N2H"'" emission; in 2A the star has already heated 
the material and some CO has been released, reducing 
the N2H+ abundance and emission; the evolved core 2B 
has been thoroughly heated by the star, and most N2H"'' 
has been effectively destroyed by the released CO. This 
hints at triggered star formation with the sources lining 
up from southeast to northwest in evolutionary order, with 
2B older than 2A, and 2A older than 2C. 

Several molecules only show emission in a very narrow 
velocity range around the systemic velocity of the cloud 
of ~ S'km s-i: T1^^C0+, SO, CH3OH, N2H+, and C^^S. 
Others show pronounced gradients in a north-south di- 
rection within « 20" from IRAS2A (HCN, HCO+, CS) 
and along the east-west outflow (most clearly in HCO+). 
Whether the north-south gradient around 2A is rotation 
in a circumstellar envelope or is related to the north-south 
outflow seen on larger scales is addressed in Sectional To 
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Fig. 13. The contrast between the N2H+ and SCUBA 
emission: the N2H+emission divided by the 450 /im con- 
tinuum emission (normalized) along a straight line with a 
position angle of 45° through IRAS2A (solid line) and the 
450 fim continuum emission in the same positions (dashed 
line). 



asses the amount of recovered flux as function of veloc- 
ity. Fig. 1141 compares single-dish spectra with interferom- 
eter spectra averaged over the single-dish beam size and 
converted to the antenna-temperature intensity scale. The 
interferometer recovers at most 17% of the emission, and 
much less in many cases. Apart from a scaling factor, the 
line shapes of C'^^S, N2H+, and H^^CO+ are similar in 
the interferometer and single-dish spectra, implying that 
although the interferometers picked up only the more com- 
pact emitting structures, only small velocity gradients can 
be present within the envelope itself. Deep self-absorption 
apparent in HCO"*", HCN, and CS near the systemic veloc- 
ity indicates that surrounding cloud material is optically 
thick and entirely resolved out. The velocity structure seen 
in these lines therefore only reflects material at relatively 
extreme red- and blue-shifts. 

4.2. Envelope contributions to the line emission 

The available envelope model can help explain the in- 
terferometric line data, by separating the expected emis- 
sion of the envelope from other components. We adopt 
the power-law envelope model of Sect. |3 and fix the 
molecular abundances from fitting the integrated inten- 
sities of single-dish observations of optica llv thin isotopes. 
Table 0] lists the derived abundances. Ijorgensen et alJ 
( 2003bi) discuss this method in greater detail and com- 
pare the results for a larger sample. With these abun- 
danc es, the molecular excitation is solved using the code of 
■Hogerheiide fc van der Ta^ l)2000.) and the emergent sky- 
brightness distribution calculated. For comparison with 
the interferometer data the sky-brightness distribution 
was sampled at the (u, v) positions of the data. The result- 
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Fig. 10. Integrated line emission from the BIMA observations for a) HCN, b) HCO+, c) N2H+and d) C^'^S plotted 
over the 3 mm continuum maps (grey-scale) . The outflow axes have been marked with straight lines with the red part 
being solid and blue part being dashed. For HCN and HCO+ the emission has been integrated over the red and blue 
parts of the line (3 to 7 km s~^ and 9 to 13 km s~^) shown as the dashed and solid lines, respectively. For N2H"'" and 
C^'*S the total integrated emission is presented, in the case of N2H+ integrated over the main group of hyperfine lines. 
For C'^^S the contours are presented in steps of 3(7, for HCN and HCO+ in steps of 5(j and for N2H+ in steps of 7<j. 



ing visibilities were then inverted, cleaned, and restored in 
a similar way as the data. 

Fig. 1151 compares the observed and modeled line emis- 
sion for C'^'^S. The upper panel shows a comparison be- 
tween the model and the single-dish observations (upper 
spectrum) and the interferometry data convolved with the 
single-dish beam (lower spectrum). In the lower panel the 
visibilities are plotted as a function of projected baseline 
length. Both comparisons show that the model works very 



well in describing the interferometry and single-dish line 
observations simultaneously and reproduces the emission 
distribution at the observed scales. This has two implica- 
tions. First, that optically thin species such as C'^^S trace 
material in the envelope and are well described by the 
model derived from the continuum observations. Second, 
for species such as C^^'S the chemistry is homogeneous 
at the observed radial scales in the envelope, so that a 
constant fractional abundance is sufflcient to describe the 
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Fig. 11. Integrated line emission from the OVRO data, showing a) CS, b) }1^^C0+, c) SO and d) CH3OH. CS is 
integrated over blue (5 to 9 km s^^; dashed contours) and red (9 to 13 km s~^; solid contours) parts of the line with 
contours in steps of 3a. For the other molecules lines have been integrated from 7 to 11 km s^^ and contours are given 
in steps of 2a. The grey-scale indicate the 3 mm continuum maps. 



Table 4. Abundances in the ID static envelope model 
der ived from of single-dish line observations. For details 
see I.Targensen et al.l ( 2003hl) . 



Molecule 




Abundance ([X/H2]) 






1.4x10"'" 


CS 




1.3x10"^ 


H"CO+ 


(all except 1-0 line) 


4.3x10"" 


H"CO+ 


(1-0 line) 


8.0x10"" 



chemistry within the assumptions and uncertainties in the 
modeling. 

For H^'^CO+ the situation is a bit more complex. The 
modeling of the single-dish lines re veal a pi cture similar 
to that of the CO isotopic species in lPaoerJ : while a con- 
stant fractional abundance can describe the line intensities 
of higher J lines, the intensi ty of the J = 1—0 line is under- 
estimated by the model. In lPaoerl it was suggested that 
this was due to ambient cloud material being picked up by 
the larger single-dish beam. The same problem may be an 
issue for Hi3C0+. Fitting the Hi3C0+ J = 1 - single- 
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Fig. 12. Fh-st order moment (velocity) maps of the a) HCN, b) HCO+ (BIMA) and c) CS emission (OVRO). Each 
map has been overplotted with the total integrated emission in steps of 3 cr (solid line contours). The outflow axes 
have been marked by lines and IRAS2A and IRAS2B continuum sources by stars. 
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Fig. 14. Comparison between the single-dish observations 
(dark) and corresponding spectra from the interferometer 
observations restored with the single-dish beam (red). The 
spectra from the interferometry observations have been 
scaled by the factors indicated in the upper right corner 
(factors 5-12) to include all spectra in the same plots. 



ferent. We emphasize, however, that this is a significantly 
smaller effect than what is seen for, e.g., CS, HCN and 
HCO"*" (Fig. I14|l . It is also seen that the model cannot 
describe the emission on smaller scales when plotting the 
visibilities vs. projected baseline length, as illustrated in 
the lower panel of Fig. ^] If the lower H^'^CO^ abundance 
found from fitting the higher J lines is adopted, the model 
perfectly reproduces the observed H^'^CO"'" emission dis- 
tribution. This discrepancy suggests that the single-dish 
beam of 44" picks up the ambient cloud and the envelope 
around IRAS2B as it also is seen in the II^^CO+ 1-0 in- 
terferometry maps. This will contribute to the spectrum 
extracted from the interferometry cube when convolved 
with the single-dish beam (e.g. the upper panel of Fig. [TB|) . 
The good fit to the visibility curve in the lower panel of 
Fig. ^1 indicates that the abundance constrained by the 
higher excitation single-dish line observations of H^'^CO'^ 
is representative of the actual envelope abundance. 

It is not possible to account for the observed CS, HCN 
and HCO+ emission within the envelope models. As can 
be seen in Fig. El the CS line intensity is, for exam- 
ple reproduced only at intermediate baselines where also 
the single-dish line observations are sensitive. On longer 
baselines the model clearly breaks down and underesti- 
mates the observed emission. It is also evident that the 
pronounced double peak structure seen in interferometry 
spectra cannot be explained with a simple collapse model 
alone. This problem will be further explored in the next 
section. 



dish line alone gives an abundance of 8. Ox 10"". With this 
abundance, the model can describe the intensity of a spec- 
trum convolved with a beam equivalent to the single-dish 
observations as illustrated in the upper panel of Fig. 1161 
The model, however, cannot fit the line profiles simultane- 
ously for the single-dish and interferometry spectra when 
a constant turbulent broadening is adopted. This is likely 
caused by the larger single-dish beam picking up the more 
extended cloud where the velocity distribution may be dif- 



5. Velocity structure beyond the envelope 

Fig. ll8lshows position- velocity diagrams of HCO+, HCN, 
and CS along the north-south velocity gradient apparent 
in the material within « 20" from IRAS2A as seen in 
Fig. I10I12I The previous section found that the velocity 
structure in these lines could not be explained by the en- 
velope model with infall. Fitting a linear gradient to the 
velocity centroid at each offset yields values as given in 
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Fig. 18. Position-velocity diagrams for a) HCN, b) HCO"'" and c) CS. The solid line indicates a linear gradient fitted 
to the centroids for the velocity channel. The hyperfine splitting of HCN is seen as the extension of emission along the 
velocity axis. 



Table 5. Velocity gradient derived through a linear fit to 
the centroid of the brightness for each velocity channel as 
indicated in Fig. El 



Molecule 


Velocity gradient 

[km s"""^ arcsec^""^] 


CS 


1.09 ±0.29 


HCN 


1.00 ±0.49 


HCO+ 


1.25 ±0.53 



Table [S] for the three species. It is seen that the fitted ve- 
locity gradients agree well and correspond to a weighted 
average of 1.10±0.23 km s~^ arcsec"^. This gradient may 
be an overestimate of the actual gradient, because the 
interferometer predominantly recovers extreme velocities 
due to resolving out. On the other hand, velocity gradi- 
ents inferred from single dish observations only are biased 
toward velocities closer to the rest velocity of the cloud 
due to the larger beam. 

The inferred north-south velocity gradient of 
1.10 km s~^ arcsec""'^, or 1.03x10"^ km s~^ pc~^, is 
two orders of magnitude la r ger th an that inferred by 
IWard -Thompson fc BiicklevI l(200ll) from single-dish 
CS and HCO"'" observations. This increase in velocity 
gradient from the spatial scales of the single-dish to the 
interferometer data is too large to be explained by differ- 
ential rotation in either a Keplerian structure (expected 
increase of a factor ^/TO) or a magnetically braked core 
ijBasiJ Il998l) (expected increase of a factor 10). Even if 
the gradients on single-dish and interferometer scales are 
unrelated, Keplerian rotation cannot explain our observed 
velocities since it requires an unrealistically large central 
mass of 31 Mq. 

An alternative explanation for the north-south ve- 
locity gradient around IRAS2A is that it is part of 
the i iorth-south outflow detected on larger scales in 
CO (iLiseau et all Il988t jEiigaj giolj^fc_^lBint^ ^^991 
iKnee fc SandelJ l2000|) . lEngargiola fc Plambeckl l)l999|) 
concluded that the origin of this flo w lies within a few arc- 
sec from 2 A, while iHodapp fc Ladd (,1995.) infer a north- 



south jet that passes within a fe w arcsec from 2 A fr om H2 
images. Neither this Paper nor iLoonev et al.l l)200fll) find 
evidence for continuum emission from a second source, 
although it could be below the detection limit or be un- 
resolved (< 0.3" — 65 AU). The different levels to which 
HCO+, HCN, and CS trace the east- west and north-south 
flows may reflect differences in shock chemistry as the 
flows progress through the inhomogeneous cloud environ- 
ment of IRAS2. This also serves as caution in interpreting 
differences in the spatial extent of, e.g., CO line wing emis- 
sion as differences in 'dynamic time scales'; this presup- 
poses similar environments in which the flows propagate. 

Instead of two perpendicular outflows, a single, wide- 
angle (90°), northwest-southeast flow could also explain 
the observed gradients. In this interpretation, what appear 
to be two independent flows actually trace the interac- 
tion of the wide-angle flow with the surrounding material 
along the sides of the cavity. This sc enario is reminiscent 
of the wide-angle outflow of B5-IRS1 ijVelusamv fc LangeJ 
.1998.') . However, in this scenario the jet-li ke morphol- 
ogy of the shocked region east of IRAS2 ijBlakel Il993 
'Bachiller et alJll998l:ljOrgensen et alJl2003a^ is difficult to 
explain. 

6. Conclusion 

This Paper has shown that the envelope model derived 
from submillimeter continuum imaging with SCUBA pro- 
vides a useful framework to interpret interferometric mea- 
surements of continuum and line emission. It allows sepa- 
ration of small-scale structures associated with the enve- 
lope from small-scale structure in additional elements of 
the protostellar environment, such as outflows and disks. 
Our main flndings are as follows. 

1. Compact 3 mm continuum emission is associated with 
the two protostellar sources NGC 1333-IRAS2A and 
2B; the starless core 2C is not detected, indicating it 
lacks sufficient central concentration. 

2. The 3 mm continuum emission around 2 A in the in- 
terferometer data is consistent with the extrapolation 
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Fig. 15. Upper panel: Comparison between the C"^ S 
emission from the single-dish observations using the IRAM 
30 m telescope with the BIMA interferometry observations 
(lower spectrum), offset along the Tmb axis by -0.8 K and 
restored with a 25" beam similar to the single-dish data. 
The prediction from the ID static m odel of the emission 
brightness distribution from lPaoer j with abundances de- 
rived from single-dish line observations given in Table 21 
and sampled at the relevant (u, v) grid has been overplot- 
ted on the spectrum in red. Lower panel: visibility ampli- 
tude plotted as function of projected baseline length. The 
predictions from the envelope model with C'^'^S abundance 
constrained by the single-dish line observations have been 
overplotted as the solid line. The zero-expectation level is 
indicated as the dotted histogram. 



of the envelope density and temperature distribution 
to small scales. Changes in the extent of the enve- 
lope and inclusion of the interstellar radiation field do 
not change this conclusion. A densi ty structur e as pre- 
dicted from an inside-out collapse (IShlJll977^ fits the 
data equally well. 
3. The 3 mm continuum data show the presence of a 
22 mjy unresolved source, presumably a circumstel- 
lar disk of total mass > 0.3 Mq. 
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Fig. 16. As in Fig. El but for H^'^CO"'' emission as traced 
by single-dish observations from the Onsala 20m telescope 
(44" beam) and interferometry with the OVRO millime- 
ter array. The H^'^CO+ interferometer spectra have been 
scaled by a factor 3 in order to be able to visualize them 
together with the single-dish observations. The solid line 
in both panels indicate the model with H^^CO+ abun- 
dance constrained by the high J line observations. The 
dashed line indicate the model constrained by the H^'^CO+ 
J = 1 — line. In the lower panel models both with a tur- 
bulent broadening of 0.5 and 1.5 km s~^ have been shown 
for the abundance constrained by the 1 — line. The in- 
tegrated intensities are the same, however. 



4. Line emission in the optically thin tracers H^'^CO"'' and 
C^'^S is also consistent with the extrapolated envelope 
model. Since these tracers are optically thin, this sug- 
gest that the bulk of the material is well described by 
the envelope model. 

5. Optically thick emission lines of CS, HCO+, and HCN 
only trace a small fraction of the material at veloci- 
ties red- and blue-shifted by several km s~^. Emission 
closer to systemic is obscured by resolved-out large- 
scale material. The detected emission is closely asso- 
ciated with two perpendicular outflows directed east- 
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Fig. 17. Comparison between interferometry observations 
and envelope model for CS. 



west and north-south. This suggests that the source 
2A is an unresolved (< 65 AU) binary. 
6. The morphology of the line emission in the maps shows 
that chemical effects are present. An example is the 
emission of N2H+ that traces cold material around 2 A 
and that is especially strong toward the starless core 
2C. The emission avoids the region around 2B and the 
outflows. We suggest that the dearth of N2H+ emis- 
sion is due to destruction through reaction with CO 
released from ice mantles in warmed-up regions. This 
indicates an evolutionary ordering 2C-2A-2B, in order 
of increasing thermal processing of the material. 

This work suggests that successful interpretation of the 
small-scale structure aroimd embedded protostars requires 
a solid framework for the structure of the surrounding 
envelope on larger scales. In this framework one can ef- 
fectively fill in the larger-scale emission that is resolved 
out by interferometer observations. The submillimeter- 
continuum imaging by instruments like SCUBA has 
proved particularly powerful because it does not suffer 
from chemical effects that make line emission measure- 
ments so complex. On the other hand, this very chemistry 
reflects which physical processes are occurring: e.g., the 
N2H+ emission that shows the thermal history of the ma- 
terial. 

The success of the envelope model in describing the 
optically thin species, such as C'^^S and H^^CO"*" makes 
IRAS2 a promising candidate in order to study the rela- 
tion between the envelope chemistry and the spatial dis- 
tribution of molecular species. In particular, studies of a 
larger sample of optically thin molecular lines at arcsecond 
scale resolution may probe difFcrcnces in the radial distri- 
butions of molecules reflecting the chemistry. IRAS2A is 
for this purpose a promising target due to the relative 
simplicity of the central envelope component. High angu- 
lar resolution, high sensitivity maps may also allow for 
a more detailed comparison to models for the protostel- 



lar collapse in order to possibly address the evolution of 
low-mass protostars in the earliest stages. 
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